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On the Inclinational Terms in the Moon's Coordinates. 

By P. H. Cowell, B. A., Fellow of Trinity College, and Isaac Newton 
Student in the University of Cambridge, England. 



In the first volume of the American Journal of Mathematics, Dr. Hill pub- 
lished some papers in which the variation terms of the Moon's orbit were com- 
puted both literally and arithmetically with extreme accuracy. When all other 
terms but those depending on the ratio of the mean motions only are omitted, 
the remaining terms denote a periodic orbit which relatively to the Sun is a 
closed curve. This orbit has been made to play the part taken by the ellipse 
in the older lunar theories. In a subsequent paper published in the Acta Mathe- 
matica, vol. VIII, Dr. Hill obtained an expression for the motion of the Moon's 
perigee, so far as this depends on the ratio of the mean motions. His papers 
were followed by others by Prof. B. W. Brown, published in the American 
Journal, in which the terms depending on the solar parallax and the lunar 
eccentricity were computed. In all these papers Dr. Hill's orbit was treated as 
" the intermediary orbit," about which the actual orbit oscillates. 

It is the object of this paper to take into account, according to Dr. Hill's 
method, the inclination of the Moon's orbit, considering it as being the mani- 
festation of a small oscillation about Dr. Hill's distorted circular orbit, which 
relatively to the Sun is a closed curve. 

The terms multiplied by the first power of the inclination have been calcu- 
lated to the sixth order, and an expression for the part of the motion of the 
Moon's node, that depends upon the ratio of the mean motions only, has been 
found as far as the eighth order; that is to say, one term further than in 
Delaunay's series. The numerical values having been found to twelve places of 
decimals, we have the means of evaluating the residue of Delaunay's series. 

The terms multiplied by the square of the inclination have been calculated 
14 
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to the fifth order ; that is to say, to the seventh order when account is taken of 
the factor y z . 

The terms multiplied by the third power of the inclination have been calcu- 
lated to the fourth order in tn . From these terms the part of the motion of the 
Moon's node that depends on the factor y % has been calculated to the sixth order 
of rn, one term further than Delaunay's series. 

The principal papers dealing with Dr. Hill's Lunar Theory are tabulated 
in a paper by Prof. E. W. Brown, American Journal of Mathematics, vol. XVII, 
p. 320. The notation of the above paper will be adopted here. 

Preliminary Modifications of the Problem to be Solved. 

Delaunay, in his Lunar Theory, introduced two modifications. He treated 
the Moon's mass as negligible in comparison with the Earth's, and the Earth's 
mass as negligible in comparison with the Sun's. It is well known that the first 
modification does not affect those terms in the final expressions for the Moon's 
coordinates which are known as non-parallactic ; that is to say, which do not 
vanish when the Sun's parallax is neglected. The terms containing the first 
power of the solar parallax as a factor can be calculated according to the 
first modification and then multiplied by the ratio of the difference to the sum 
of the masses of the Earth and Moon. The effects of the first modification in 
the terms containing the second and higher orders of the Sun's parallax and the 
effects of the second modification — treating the Earth's mass as negligible in 
comparison with the Sun's — cannot be allowed for so readily, but it may be 
inferred that, if the solar parallax be always affected with the factor denoting the 
ratio of the difference to the sum of the masses of the Earth and Moon, no sen- 
sible error is introduced by the first modification, while the second modification 
may be adopted and its results accepted as they stand without sensible error. 

For these reasons the lunar theorists, with the single exception of Hansen, 
have modified the problem in the above two ways before entering upon the 
analysis. 

The Differential liquations. 

Taking the axis of a perpendicular to the plane of the ecliptic, and the axes 
x, y moving in the plane of the ecliptic, so that the axis of x is constantly directed 
toward the mean place of the Sun ; also taking as a special unit the time in 
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which the mean elongation of the Moon from the Sun increases by one radian ; 
then the effective forces for the Moon are 

d?x n dy o 

W~ 2m df- mx ' 

J?]L + 2m^—m % v 
df + dt my ' 

d?z 
df ' 

where x, y, z are the coordinates of the Moon relatively to the Earth, and rn is 
the Sun's mean motion, or, which is the same thing owing to the units of time 
adopted, the ratio of the Sun's mean motion to the difference of the mean 
motions of the Moon and Sun. 

Denoting by x the mass of the Earth, and in accordance with the second 
modification, putting wiV 3 for the mass of the Sun alone, instead of for the sum 
of the masses of the Sun and Earth, where a' is the Sun's mean distance, the 
impressed forces are derived from a potential function £l where 

x , m z a'' mV o 

~ V(r 2 + s 2 ) + V(r'*— 28/ + r a + s a ) "T 7 ^ ' 

where Prof. E. W. Brown's notation, reproduced below from Amer. Jour., vol. 
XVII, p. 321, is adopted : 

2a', r', e' the major axis, radius vector and eccentricity of the Sun's orbit, sup- 
posed elliptic ; 
x, y, z the coordinates of the Moon ; 
v the solar equation of the centre ; 
£= x cos v + y sin v ; 
n'^a!* the mass of the Sun ; 

x the sum of the masses of the Earth and Moon ; 
r* = x 2 + y % . 

We shall subsequently introduce 

u = a; + yV — 1 , s = x — yV — 1 , 

i^-e , »-<>$£ ~ s/=\ dt * 
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Forming the equations of motion, and taking the centrifugal terms over to 
the right-hand side, we have 

x-2my = ^-, 

y + 2mx — -g— , 

— ar 

where 

O! — £l + \m % (a? + y 2 ) = £1 + £mV. 

Transforming to the complex variables just defined, the equations become 

dCl' 



B*u + 2rnDu= — 2 
IPs — 2mDs — - 

Dh =- 



ds ' 

da 1 

du ' 
d£l< 

dz ' 



where 



m z us 



£1' = , * ,.. + mV 8 1(/ 2 - 2Sr>+ us + z*)'* - S.r>- 3 \ + h 

X ,12 i 3 «'* 3S 2 —US Z 2 

(«* + z 2 )* r' s 2 

, m 2 a' 4 5^— 3S(us+f) , 

+ W T 5- 2 +••••, 

omitting the term —j- which disappears on differentiation. 

From CI' we select those terms that remain when the solar eccentricity (and 
therefore the solar equation of the centre) and also the solar parallax are put 
zero. In this case #= x = (u + s)/2, and the terms selected are 

-. — * + f «i 2 (« + s) 2 — imW. 
(its + sr) 1 

Following Prof. E. W. Brown, we put 

fl' = « + 1 m 2 (« + s) 2 — *mV+ *"i , 
(«s + 2r)' 
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where 



fl, = 3m 2 [j£ & - *(« + s) 2 ] — m 2 (tM + s 2 )^ -l) 

+ .?£. . ?lL [5S 3 — 3S(us + z*)] + .... 

The terms which will be investigated in this paper are those that arise 
from treating the third coordinate of the Moon as small in comparison with the 
other two, instead of entirely neglecting it, as Dr. Hill did, in investigating his 
variation curve. The solar parallax and eccentricity will be neglected, and 
consequently fi x is neglected also. The differential equations therefore become 

(Z> 2 + 2mD + |m 2 ) u + f m*s - . *" = 0, 

(us + zy 

(P- 2mD +|m 2 ) S + |m 2 «-^i-^ = 0, 

(US + 2 2 ) S 

When z is neglected, the solution involves four arbitrary constants : of these, 
two, viz. the linear dimensions of the orbit* and the time of conjunction, occur 
in Dr. Hill's investigations of the variation curve. This curve represents a 
periodic motion, but inasmuch as two arbitrary constants are still wanting, it 
is not the most general solution of the equations in their modified form. The 
general solution should be regarded as denoting an oscillation superimposed upon 
Dr. Hill's periodic motion. The constants introduced define the amplitude of 
this oscillation and its phase. They replace the lunar eccentricity and the time 
of passage through perigee of the theories that take as the intermediary the 
ellipse and not the variation curve. As however in this paper we do not pro- 
pose to retain the lunar eccentricity, we quote from Dr. Hill's papers (Amer. 
Jour., vol. 1) the particular solution of the differential equations 

z— 0, 
« = w = a 2a i £ 2i + 1 , 

* — *o ~~ z a o^ a — *— is i 

*m may have any arithmetic value consistent with convergence. Hence it appears that the mean 
motion and the mean distance are both arbitrary. This is because k, the mass of the Earth and Moon, 
will be eliminated before integration. It follows that « is a definite function of m and a<, (the linear 
dimension), and possibly of other constants as well. 
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the algebraical and numerical values of the quantities a t being given in the 
papers referred to. 

Let y denote a constant of the same order as z . Then we are seeking the 
most general solution of the form 

Z ^— Zy "f~ Zyi "P ZyQ ™f" • • • • , 
U = U Q + Uy, -\~ Uy4 -^ • • • • , 
S ^^ Sft -J" SyZ "J" Syi "J" • • • • J 

where the suffixes denote the order of the terms in y. 

The most general solution should contain six arbitrary constants. Two, 
however, disappear, as we are not considering the lunar eccentricity. We 
should therefore still have four. Besides the two introduced by Dr. Hill, we 
shall have one, K, replacing Delaunay's y, defining the amplitude of the oscil- 
lation about the variation curve and another defining the phase of that oscilla- 
tion at a given epoch. 

It is useful to regard the terms involving K as defining an oscillation about 
a periodic motion, for we are thereby enabled to predict the character of the 
solution. So long as the oscillation is small, so that the square of its amplitude 
may be neglected, we know that its period is independent of its amplitude. The 
principal part of the motion of the Moon's node is therefore a function of m only. 
Considerations of symmetry show that u, s will involve only even powers and z 
only odd powers of K in accordance with the type of solution assumed above. 
The period of the oscillation, more commonly called the period of revolution with 
respect to the node, will clearly only involve even powers of K. 

The First Power of the Inclination. 

We must try and satisfy the differential equations accurately to order K 
with the assumed solutions 

u = u , s=s , z = z y . 

The first two equations are obviously satisfied. The third equation, on neglect- 
ing K 3 , becomes 
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The solution is of the form 



Zy ■■ 



where K = 1 . 

The quantities Kj are all determinate : K is one of the two arbitrary con- 
stants required : the other has not been put in evidence, but corresponding to 
the definition £= e t/ ~ 1< - t ~ t °\ we may suppose £ ? = e v ~ lg{t ~ ti \ in which case t x is 
the second arbitrary constant required, defining the time of passage through the 
node. No confusion will be introduced by these definitions, since g is incom- 
mensurable, and integral powers only of £ and t/> are required. 

The expression m 3 -j- % ., involves only quantities calculated by Dr. Hill, 

and which we therefore suppose known. 

Following Prof. E. W. Brown's notation, we write 

The quantities M t are calculated arithmetically by the method of special values, 
and algebraically by a process explained by Dr. Hill (Acta Math., vol. VIII, 
p. 1 2). Dr. Hill has not however obtained the algebraical values to the order 
required for this paper. By his methods the following results have been 
obtained : 



5 , . „ , 9 , , . „ , 34 .,.,„,, 2704801 



2M =1 + 2m + -^--m z + 0.m 3 — -^-m 4 + 4m 5 -f ^-m t +15m 1 + 



2 ' ' 2 5 3 2 13 . 3 



13 o3 



m° 



„,, 3 „ . 19 3 . 20 4 . 43 5 , 18709 . . 759413 . 

2M t = -2~ m + -W m + T m + W + ~2KW + 2^¥^ ' 

„,, 33 4 . 2937 5 . 23051 6 

2^= -F m+ ¥X m+ 2l3^ m ' 

2Jf 3 = i^««. 

Generally M t is of order m zi . 

Substituting the assumed solution in the differential equation, we obtain 



(2j + gfK j -2^K i M j _ t =0. 
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By giving j the values 0, ± 1, =fc 2 ±k, we obtain 2& + 1 equations 

which will determine the ratios of K s to K Q (/= ± 1 , . . . . ± h) and also the 
value of g. In the equation above written the coefficient of K$ is (%+gf — 2if 
in general, a large quantity, whereas every other coefficient is of the second 
order at least. The equation therefore serves primarily to determine K; . If 
we put j = we get 

g 2 =2^M_ i K i jK , 

i = —oo 

an equation for determining g. 

From the last equation we obtain 

g* = 2M to order m 2 , 

whence , , . 3 « 

g = 1 + m + -y m 2 . 

Hence, whereas (2/ + </) 8 — 2M is generally of order m°, when / = — 1 it is of 
order m. The equation for K_ lf 

[(- 2 + gf - 2M ] K_ x = 2^K t M_ x _ u 

— oo 

cfc-i 

shows that E_ l .K is of order «i. The general law that K_ k :K is of order 
m 2 * -1 , and that .fi^ : K is of order*m 2fc now follows easily. 
Proceeding by continued approximation, we obtain 

ii . 3 a 33 3 105 4 . 43 5 , 2567 . 



, 347699 , , 6442309 . 
+ 2 16 .3 2 + 2 18 . 3 3 ' 



and the subsidiary quantities 



2 i i „ , 5 , 9 3 201 4 3221 , 

, 1031 6 . 917555 , , 13397317 8 
+ ^3 ffl+ ~2^F- m+ 2»3 3 W ' 
. n , f 9 3 183 4 7317 5 138233 ,. 

3506125 7 _ 9959833 
2 15 .3 2 m 2".3 2 ' 
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(2 + gf— 2jf = 8 + 4m + 3m z —^-m 3 — ~m i , 



3 



(— * + gf— 2if = S — Sm—6af + ±?*m l , 

(4 + gr) 8 — 2Jf = 24 + 8m + 6m 2 , 

(— 6 + gf — 2M = 24— 12m, 

(6 + <?) 2 -2Jf = 48, 

and also the coefficients 

tt tt 3 29 a 2029 . 18875 , 44633 . . 5569883 - 

2 2 2 3 2 3 2 3 2 3 

e>- j?- 3 o . 1 3 , 197 , . 3067 5 348449 6 

K i: K = _ m2+ _ m s + ^m* + ¥ r ¥ m^ - -^^ ««, 

k- k- 9 o 105 4 8553 5 . 123367 6 

K_ 2 :K = -$r m --W m ~WJ + 2 15 .5 2 m ' 

tt ir 25 4 « 803 5 . 377701 , 



2 8 ^ 2 7 .3.5 ' 2 U .3 2 .5 2 

75 5 18797 

3 * J1 "°~ W m ~ 2 1S X5 

833 
2 12 .3 

Comparison with Delaunay's Expressions. 



-tr xr 75 5 18797 « 

K_ S :K = — — m 5 — gl3 - - m 6 



From Dr. Hill's expressions for 'a; and y we form (a; 2 +# 2 )/ao, and hence 
a (a 2 + s/ 2 ) - *. Multiplying by z/a , we get an expression for the tangent of the 
latitude, or correctly to the first power of the inclination the latitude itself. 
This expression for the latitude is 
n „rA 3 s 141 4 3893 . 77345 A . „ 

, /ll a . 13 3 . 943 4 . 6017 . 3148129 6 \ • ,„ n , m 

. /161 4 . 3607 5 . 1282109 A • /An , ™ 

15 
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. / 3 13 2 1133 3 13243 4 

V 2 s 2 6 2 9 .3 2 n .3 a 

17617 5 . 6400195 A . , on , ™ 

~ ~¥K¥ + 2 16 ,3* m J sin ( ~ + } 

, / 33 3 225 4 50563 . 75041 A • i „ n . m 



. / 483 5 83333 A . , .„ . „,-] 



This is easily seen to be identical with Delaunay's expression (Delaunay's rn 

denotes the ratio of the mean motions), provided that we equate Delaunay's y 

^ vTi 3 s 141 4 3893 . 77345 {] ™. ,, . 

to K II — —j nr -g- m 4 — -jg— - m 5 — 18 % m b . This result is necessary 

to enable us to compare the two theories. 

Alternative Calculation of g. 

The quantity g may be determined alone without simultaneously calculating 
the quantities K. The method pursued is analogous to that in Dr. Hill's paper 
on the motion of the perigee (Acta Math., vol. VIII). 

If we eliminate the quantities K from the equations 

l(2j + gf-2M ]K j -2^K { M j _ i = 

1 = 00 

we obtain a determinant of an infinite number of rows and columns. 





1 


—2M 1 


2-Jf a 


2ifcT 3 


2Jf 4 


2M t 


(i+gy-2M 

1 

2M X 


(4+0) 2 -2M o 
2Jfj 


(4+g)*-2M 

2Jf 2 


(4+sr) 2 -2if •• 

2Jf 3 


' ' (2+g)*-2M ' 
23f 2 


(2+flf) 2 -2Jf ' 
1 
2Mi 


(2+fl-) 2 -2ilfo ' 
2Jfj 


(2+90 2 -2Jlf o • • 

2M 2 


g 2 -2M ' 
2M 3 


g 2 -2M ' 

2K 2 


g 2 -2M ' 
1 
2Kj 


g*-2M 

2M 1 


' ' (— 2+g)*-2M„ ' 
2M 4 


(— 2+sr) 2 — 3ifcT ' 

2M 3 


(-2+gy-2M ' 
2ilf 2 


(-2+gy-2M •• 
1 


" ' (-4+0) 2 -2.af o ' 


(_4 +g ,) 2 _2i)f ' 


(-4+g) 2 -2M ' 


(-4+3) 2 -2K ' 





This equation determines g. Its use, according to M. Poincar6, is legitimate, 
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since the sum of the elements not situated on the leading diagonal is finite, and 
the determinant therefore converges absolutely and uniformly. 
One solution of the equation 



Dh: 



(f +m2 > 



is z = 2K^+°. 

This solution may be written 



Z : 

i 



when Kf = K i+j . 

In this second form of the solution, g -\- 2i replaces g, and consequently we con- 
clude that if g is a root of the determinantal equation, g + 2i is also a root when 
i is any integer, positive or negative. This result also follows from an inspection 
of the determinant. Since also the determinantal equation is unaltered on 
changing the sign of g, we conclude that if g be one root, ±2i± g (i, integral) 
is also a root. There can be no root not included under this expression, for 
else there would be a third solution of the third differential equation of motion 
linearly independent of the two known solutions 

z = 2K j ^+» 
and z = 2 %£-»-', 

which is impossible for a linear differential equation of the second order. 

Since the roots are not repeated when we limit ourselves to a finite number 
of rows and columns, we infer that the roots are not repeated in the case con- 
sidered. The infinite determinant must therefore be identically equal to 

A (cos ng — cos 7i</ )/(cos 7ig — cos its/ 2M ) , 

the denominator representing to a numerical factor the product of the denomina- 
tors in the rows of the infinite determinant. When a finite number of rows are 
considered, the coefficient of the highest power of g is independent of the quan- 
tities M tt and we infer therefore that so is JL On putting M t = (i :£ 0) , we 
then obtain A = 1 . Finally we put g = in the identity, and we obtain 
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sm 2 — g 


D(o) = 








sin s -J-V2ifo 




i 

2M X 

" 2 3 — 2M a ' 

2M % 

O 2 — 2M Q ' 

2M B 

2»—2M 

2M % 


2Mj, 

4 2 — 2Jf ' 

1 , 

2M X 


2M i 

4*—2M ' 

2M 1 
2 i —2M ' 

1 , 

2M 1 

2 2 — 2 if ' 

2if 2 


2Jf 3 

4 s — 2M ' 

2Jf 2 

2 s — 2if ' 

2^ 

O 3 — 2ilf ' 

1 

2M 1 
2 2 — 2ilf ' 


2M i 
4 3 — 2 Jf ' " 

2M 3 

2 2 — 2if ' " 

2M % 


3 — 2M Q ' 

2M % 

2*—2M ' 

2M S 
4 2 — 2 Jf ' 


O 2 —2M ' " 

2M 1 

2* — 2M ' " 

1 , . . 


" 4?—2M ' 


2 2 — 2^' 





In this equation we may drop the suffix after g . 

Dr. Hill has developed a determinant to □ (0) (Acta Math., vol. VIII, 
p. 29-31). Using his results, but suppressing terms of order m 10 , we obtain 

n (0) - 1A * COt (^ Vm) \ ^ AMI } 

Ulj ~ 1+ W(2M ) )l — 2Jf 4-2Jf j 

ncot(^y(2M )) t(nV{2Mo)) _ 1 a ^ W l} 4 

^ S2V(2M )(1— 2M f} V(2M ) 2M ^ l-2if ^ 2.(4-2if )J v u 

37t cot (-5-^(21^)) 

+ V 2 — <~-^ (2M 1 f. 2M, 

+ 8 V(2Jf )(l - 2Jf )(4 - 2M Q ) K l} 

n cot (-|- V(2Jf )) , j^ j^ 9 \ ?tcot7tV(2ilf ) 

+ 128V(2if )(l-2if )4 V~ 2iT + 1 - 2Jf„ + 2.(4-2if )J ^(2if ) 

25 1 2_ 4 9 

_ 16Jf + AMI + 1 - 2Jf + (1 - 2Jf ) 3 8 . (4 - 2M ) 

, 9 4 7t 2 Um \6 

+ (4 — 2Jf ) 2 ~~ 9 - 2M ~ 6Mj K V ' 

The last term has been written down because some parts of it are apparently of 
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the eighth order. These parts, however, cancel out, and the term is zero to the 
ninth order inclusive. 
Now 



^-f~9 l-rin»-S-(f-l) 



2 v v u ' " 2 



sin 8 ^ V(2M ) 1 - sin 2 A (V(2if ) - l) 

= 1 -4{ ( ^- 1) ~^ 2J/( » ) - 1 ) 2 } + -" 



Hence (V(2Jf )— l) 2 — (g— If = coefficient of — - in D (0) developed in powers 

of n. 

To obtain the coefficient of 7t 2 in D (0) , 

cot -J- V{2M ) may be replaced by — -J- (\/(2Jf ) — 1), 



ncot(nV(2M )) " 



V(2Jf )— 1' 
and 

3 „ 3 . . 21 , , 143 



,5 



V (2iI/ ) = l + m + ^r ™* - -gr ™ 3 + ~Y m * + ~*~ m ' 

^V3 '"' "*" I^T "' " r 2 U .3 

Hence to the coefficient of — in D (0) , 

4 



2231 , , 2431 , . 1058251 . 
+ irnr m + -otnr m + OUQ3 m 



37t cot (^V(2MS) 297 59103 

V Zi ' (2J^) a 2il^ contributes ™- rn' + ^i^ < 



8V(2Jf )(l — 2if )(4 — 2Jf ) 2 s ' 2". 5 



ncot(^V(2M )) ^ {2Mo) 



1 + 2 + 



3 2 V ( 2Jf ) (1— 2Jf ) a L V ( 2 Jf ) 2 Jf ^ 1 — 2M r 2 (4 — 2 Jf ) -1 

, 1215 8 24165 
contributes — „ 13 m 8 ^ir - m > 

4 V(2Jf ) 4 — 2M 2 9 

ncotj^W ) 

i contributes 



4V(2if ) 1 — 2J/ 2V(2if )(V(2i/„)+l)- 
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Now 

.„, 9 4 . 57 . , 681 , , 233 , , 432661 8 , 19693 , 
AMt = ^ r rn i +^ r rn* + -gy- tn* + — m 7 + ^g, to 8 + -^-^ w 9 , 

2V(2if )(V(2M ) +1) = 4+6m+ -^ m» - A m s + -|- m 4 + ^ »»», 

whence the term 

_ 9 4 , 87 5 . 723 6 . 2587 , . 23941 8 , 42469 , 

- -¥ m + ~w m + ~r w 6 + ■ ¥T m ' + -^ « 8 + -^rg •» 9 - 

Also 

(V(2if )-l) a = m a + A m 3 _ _^_ m4 _ ^5_ ffl5 + 707 m6 

, 3329 , . 91475 8 . 1468975 , 
+ 2*73 w + 1^ TO + ^¥~ m - 



Hence 

V» = ™ 2 -i-- „,, „, ,. 

2 2 2 4 2 



(q — 1 Y = m* A m 3 m 4 — — j- m 5 r wi 6 



,371 . , 283 . . 135739 „ 

+ 2^3 W +^ F m +^^ m ' 

and hence 

1 , . 3 . 33 3 105 , . 43 . 
g = i + ffi + _ m* — -g- m % - -^ m 4 + -^ir ™ 5 

, 2567 6 . 347699 , , 6442309 . 

the same result as before. 

This agrees with, and is carried one term beyond, Delaunay's result. 
The coefficients in the series in powers of Dr. Hill's m are considerably 
smaller than the coefficients of Delaunay's series. The increase in convergency 
is even more striking in the expression for the latitude : for instance, the term 

1955563 ._ 
2 16 .3 4 
100407473 



y .»i 6 sin ( — 2D + F), which here has a coefficient — --^ — r ~ in Delaunay's series 



has a coefficient — 



2 ie 34 • 

Arithmetical Results. 



Dr. Hill has calculated by the method of special values an expression for 
%/r* with numerical coefficients (Amer. Journ., vol. I, p. 249). We hence obtain 
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2M = 


1. 


17804 


45712 


77166 


2M,= 


0. 


01261 


68462 


48930 


2M % = 


0. 


00012 


57766 


75006 


2M 3 =z 


0. 


00000 


12059 


39899 


2Mi = 


0. 


00000 


00113 


02925 


2M 5 = 


0. 


00000 


00001 


04375 


2M e = 


0. 


00000 


00000 


00954 


2M 1 = 


0. 


00000 


00000 


00008 



The following arithmetical values have been obtained : 



g=l. 08517 
id = —.03698 39313 94 
K_ % = — .00004 65750 01 
JT_ 3 = — .00000 01755 37 
K_i = — .00000 00008 8 7 
ir_ 5 = — .00000 00000 05 



14265 58 

K x — + .00151 22192 28 

K % = +.00000 58 673 61 

JT 3 = + .00000 00299 8 2 

JT 4 = +.00000 00001 75 

K 5 = +.00000 00000 01 



In performing the arithmetic, the coefficients of the quantities K were 
replaced by their arithmetical values. Since the arithmetical values of the 
quantities K were already approximately known from the series representing 
them, it was obvious to how many decimal places each coefficient had to be 
calculated. 

To this value of g corresponds a motion of the node given by 



n 



** = g- 1 -** = .00399 
at 1 + m 



91645. 



Delaunay's series gives 



.00399 91772, 



a difference equal to about .00000 3 of the whole. 

As the error between theory and observation amounts to .0006 of the 
whole, it appears that it cannot be accounted for, as Delaunay supposed, by 
the incompleteness of that part of his expression which depends on the ratio of 
mean motions only. 

Adams (M. N. R. A. S., Nov. 1877) has obtained the result : 

g — l 08517 13927 46869. 
The discrepancy between the value obtained in this paper and that obtained by 



114 Cowell : On the Inclinational Terms in the Moon's Coordinates. 

Adams is due to the different assumed values of the mean motions. Adams uses 

«l = 0. 08084 89030 52, 

whereas in this paper Dr. Hill's value 

m = 0. 08084 89338 08312 
is employed. 

It is easy to verify that the difference in the two values of g divided by the 

difference in the two values of m is the arithmetical value of dg/dm as calculated 

from the algebraic series for g. 

Second Powers of the Inclination. 

We assume as the solution of the differential equations 

u = u + u y , , s = s + s y , 

correctly to order y 2 where 

3 

and ?.«,. = aoK^[(7ck)j $-*-*'+ (*%)£-* + (##)i£~* + , \l • 

3 

The differential equations to be satisfied are 

(Z> 2 + 2m D + | m 2 ) u + fm 2 s— . * w ». = 0, 

(■P 2 - 2wiZ> + fm»)« + **»*"- (w + ^i =0. 



where zV— 1 = JEa^ (£* + * — f*-») 

and satisfies the equation 

(JP-m>)z — g =0. 

V 7 (MS + 2 8 )' 

The quantity (ws + a 2 ) -1 is inconvenient to deal with, and must be eliminated. 
First, multiplying the equations by s and u and subtracting, we get 

D (— uDs + sDu + 2mws) — f m 2 (tt 2 — s 2 ) = . (i) 
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Again multiplying the equations by s, u, 2z and adding, 

IP (us) — 2m (uDs — sDu) + f m s (u + sf — 2Du Ds 

+ 2zD\ — 2m*z* — , — ^-^ = . 
(us + a 2 )* 

The last term must be eliminated by the help of Jacobian integral of relative 
motion, obtained by multiplying the equations by Ds, Du, 2Dz, adding and 
integrating. 

Du Ds + Dz 2 + I m 2 (u + sf — <m V + -. — ?£--, = C , a const. 

(ws + zry 

Adding this equation to the last, 

D 2 (us) — 2m (uDs — sZ>«) +fm 8 (» + sf— DuDs = C + 3mV— 2aZ> 2 z — Dz\ (ii) 

Calculation of the Quantities (Htf). 

If we substitute the above assumed values in the equations (i) and (ii) and 
equate the coefficients of £** to zero, we obtain the two equations 

Y(M') r 1 4 i 2 + ( 2 *' + 1 - 2 i)( 2 *+ 1 ) + 4m ( 2 *+ 1 -i)+* m3 l^-i 

i L + {4/ + (2» + 1 + 2/)(2t + 1) + 4w(2* + 1 + /) +* m 2 \ a i+j 

+ %m 2 a j _ i _ 1 + %m 2 a_ j _ i _ 1 'J 
+ 22^+^14/ + 2y(2» +^) + (2t + gf- 3m*\ = 0, 

i 

and 

4; 2(2* + 1 — j + m) a t _ s (W) t + V^C 2 * + 1 +,/ + «) a i+j (Md) t 

i i 

— 3m»2«*-*-i (*#)« + 3«i»2 a -i-'-i^)'= °- 

In these equations (M^ and (kJd)-j occur with coefficients that are large in com- 
parison with the coefficients of the other unknown quantities. Hence the equa- 
tions may be considered as determining (Mf)j and (M')_j, and we may replace 
them by two others in which (Mc!)_j and (Mf)j have been respectively eliminated. 
The equations thus obtained will be found to differ in the sign of/ only, and one 
16 
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equation may therefore be suppressed, on the understanding that the other holds 
for negative as well as positive values of/. We are thus left with the following 
equation : 

- 4/2, U - *)[[/] - U (- 1 +j - m)] a_ s (W) s -* 
-4/^[[i]-4(/-»)(-l+/-m)]a i (^ + i 

+ ^f- [If] - 4 (- 1 +/- m)(i + 3/+ 2m)-] 2,a_ i+J _i (*#)< 

+ ^f 2 [-[i]-4(-l+/-^)(-l+3/-2m)]2 i a_ i _ j _ 1 (^) i 

-8/(-l+i- m )2{4/ a -2/(2i + 9r) + (2i + 9 r) 2 ~3m a ^^- j 
= 0, 
where [/] = 8/' — 2 — 4m + m 2 . 

This equation is an identity when/= 0. Hence the number of the quanti- 
ties (hid) exceeds by unity the number of equations to determine them. We 
may therefore assign any arbitrary value to one of them (JkTd) . Now 

is the mean value of the projection of the radius vector upon a line drawn 
towards the mean position of the Moon. a has already been supposed to be 
equivalent to this quantity, when powers of K are neglected. Any assumption 
as to (M/) is in reality equivalent to a more accurate definition of a . 

Let us therefore assume (M') = 0, an assumption that simplifies the equa- 
tions considerably, and which is equivalent to assuming that a is equal to the 
mean value of the projection of the radius vector upon a line drawn to the Moon's 
mean place, accurately, even when quantities of the order of the square of the 
inclination are taken into account. 

Note. — We might use the arbitrary condition at our disposal so that the 
relation between h, the mass of the Earth and Moon, and a and m remain the 
same as that found by Dr. Hill (Acta Math., vol. I, p. 145) even in the higher 
order of approximation which we have now reached. Otherwise the relation 
referred to must now be modified by the addition of terms proportional to the 
square of the inclination. 
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Solving the equations, we obtain 

nin 3 , 29 . , 661 3 26125 4 1385395 5 
(**0-i = -# m + ^r m + $T3 m ~ ^ir^r m 2 .3 j3 3 m > 

/77 ,, 3 , 23 3 109 4 29857 . 

(kV\ = ~^r m W W^ + ~¥K¥ ' 

W-.= *r« + f »' + S>»» 

/77/\ ■"" 4 oJiyy k 



(M0_3 = -ot »»' 



2 8 . 3.5 

75 

2 io 



correctly to order m 5 . 

The arithmetical values of these quantities are added. 

(&#)_!= + . 07279 42, 

(M/)i = — .00282 74, 

(M/)_ & = + .00009 27 , 

{W) t = — .00002 19, 

(M')_ 3 = + .00000 04, 

{kkf) 9 = — .00000 02. 

In calculating the arithmetical values, we use the equations for (&#),• and 
{hH)-j first obtained, as the labor in computing the coefficients is much less than 
for the equation when (TcH)-j has been eliminated. $H)-j can then be com- 
pletely eliminated, whereas in the algebraic equation it still has the small 
coefficient 

4i 2 l[i]-8i(-l+y- m )}a_ 2 ,+ ^|-|[y]-4(-H-y-m)(l + 3iH-2«)}a 2 ,_ 1 

+ ^f-{-[/]-4(-l+y-m)(-l + 3y-2 m )}a_ 1 , 

the last term being of order »n 4 . 

Moreover the quantities K do not appear in the second of the two simpler 
equations. 

Calculation of the Quantities (M) and (Jc'k 1 ). 

On substituting the assumed values of u and s in the differential equations, 
we obtain from the coefficients of ^ +2g the two following equations : 
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2{4(j + gf+ (2i + 1- 2j)(2i + l+ 2g) + 47n(2i + l— j+g) +lm>\a l _ j (Ick) t 
+2|4 (J + gf + (2t + 1 + 2j)(2i+l — 2g) + 4m (2t+ 1 +/— g)+ %m 2 \a i+J (h'Td), 

+ im^a_ i+j _ 1 (ltk) t + *m»2«L < _ i _ 1 (##), 
=2{(2* + </) 2 + (2/- 2* + sf + (2* + g-)(2y— 2» + gr) - 3m 2 }iT^,_ s , 

and 

4 U + g)[2(2i + l-j+g + m) at-i @k)i + S (2t + 1 +y - gr + m) a i+ ^ (#&')*] 
— 3m a 2a_i +j _i (M)< + Sm i Xa_ i ^. j _ 1 (VK) t = . 

The terms with the largest coefficients are (Jck)j, (k'h')_j. Eliminating each 
of these quantities in turn, we get the two equations 

— 4 (y + g) Z (j + g- •)[[/] - 4» (- 1 +y + g - to)] a_, (**),_« 
-4(y+ 9 r)2i[[y]-4(y+ 5 r-i)(-l+y + 9 — ro)] a, (##)-.» + < 

+ ^p [if] - 4(- 1 +y + gr-TO)(l + 3y + 3gr + 2m)]2 fl L l+J _ 1 («*) < 

+ ^[-Dl - 4(- 1 +j+g-m)(-l + 3y + 3gr-2TO)] Sa_«_, _!(##)-, + « 

+ 4(y + ^)(-i+y + ^-TO)2{(y-2*) 3 + 3(y + ^~3TO a f^^_ i = o, 

and 

- 4 (y + g) X (j + g — iW] - 44(1 +y + g + to)] a { (##)_,+, 

— 4(y+ sr)2t[[y] - 4(y+ flr-»)(i +y + fif + «*)] a-.^ 

+ ! w 2 [[y] _ 4 (1 +y + gr + to)(- l + sy + Sflr — 2m)] Sa^.^ (##), 
+ f ™ 2 [— [y] - 4 (1 +y + <? + w)(l + 3y + 3«7 + 2m)] 2a_ i+J -_ 1 (&&)« 
+ 4(y + 9 r)(l+y + gr + m)2|(y-2^ + 3(y + 9 r) 2 — 3m 2 } ^^,= 0, 

where [y] = 8 (j + #) 2 — 2 — 4m + m 2 . 
Solving the equations, we obtain 

/77N 3 . 51 s . 543 4 . 2277 5 



/7,7A -, , n 19 2 , 209 3 5109 4 , 89363 , 



2 3 2 5 2 9 2 

19 , . 209 3 5109 t . 89 
¥ rn> + ^ r rn°-^m*+- 1 ^ 

(^)i = -W m 2V3T5 m ' 

(^)-i= \rrf + -| m* + J^ m* + J^L. ««, 

59 , 2451 3 . 41541 4 2595 
^r™ 2 2r-^ 3 +-2ir-^ 4 gff 



/77X 3 . 159 , 2451 3 . 41541 4 25955 . 
(M)_ x = - — to 4- _g- to 2 -g- m 3 + — -jj- to 4 -jj— to 5 
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(##), = -^m- 


83 o 67 q 
2 5 ^ 2 6 .3 


228125 4 2450305 „ 
2 U .3* 2 13 .3 3 


(##)-.= 




25 4 . 803 B 

— =- mr A 5 mr, 

2 8 ^ 2 7 .3.5 ' 


(M)_ 2 = 




71 4 861 5 
2 n m - 2 u "» 5 . 


.(*#)• = 


27 3 

2 6 


93 4 13 6 

2 8 2 9 .5 ' 


(M)_ 3 = 




27 4 , 171 5 - 
2 io m + 2 n m ' 


(##) 8 = 




375 5 

2 io 


correctly to order m 5 . 






The arithmetical values are as follows : 






(M) = + .00165 


68, 




(M)i = + .00001 


41, 




(hh\ = + .00000 


01, 




(M)_! = — .09302 


24, 




(M)_ 2 = + .00080 


37, 




(M)_ 3 = + .00000 


14, 




(##)„ = + .98752 


58, 




(##)_!= + .00150 


88, 




(W)_2=+ .00000 


59, 




(##)■! = + .04328 


71, 




(##), = + .00020 


85, 




(##), = + .00000 


14. 



A remark similar to that in the case of the quantities (M!) may be made about 
the arithmetical work. 

Comparison with Delaunay's Results. 

Prom the values of u, s an expression for the longitude may be calculated. 
On comparison with Delaunay's expression (where, however, the symbol m has 
a different significance), we find that the two expressions are rendered identical 
if we replace Delaunay's y 2 by 



jp O-t"'-t«'-St"0 
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This relation between the constants of integration of the two theories agrees 
with the similar relation previously obtained from a comparison of the expres- 
sions for the latitude. 

Third Powers of the Inclination. 

We now return to the z equation, 

(JP-m*)z- r -Jj* =0, 

(us + g*)* 



but, since the factor (us + z a )~ f is inconvenient, we eliminate it by means of one 
of the other equations of motion, 

v ' (us + tfy 

Multiplying therefore by u, z and subtracting, 

D (uDz — zDu) — 2rnz Du —\<m\ (5u + 3s) = 0, 

and we assume as the solution 

»S—1 - a K2Kj (£*+« — $-*-«) 

+ 2(0), (£»+*— ^*-*)]. 

The motion of the node does not depend solely upon the ratio of the mean 
motions, and the above equations can only be satisfied (to order K s ) by assuming 



9 = 9o + 9^ K ' 






where g is the value of g already calculated and independent of K, and g^TP the 
correction depending on the square of the inclination. 

Calculation of the Quantities (#%% 

Substituting in the differential equation, and equating the coefficient of 
£%+i+g t zer0) we obtain 
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2[(2i + g + g y ,Ky—(2j-2i + iy-2m(2j-2i + l)-% m *]K i a j _ i 

— f m 2 22T i a i _ ; -_ 1 
■ pfS l(2i + g f—(2j- 2i+ If- 2m(2j- 2i+ 1)— f m 2 ] 25 (&#),_, 

+ 2 [- (2i + gr ) 2 + (2/+ at +1+ 2^ ) a +2m(2y + 2» +1+ 2<? ) + *«»] K t (kh) J+i 

+ §m*2K i (Jj'jc!)_ i _ j _ 1 
+ 2 [(2» + g f - (2/ - 2t + l) 2 - 2m (2/— 2t + 1) — f O (#»#), o,_, 

where 2 includes all integral values of i, positive, negative and zero. 

In this equation we may equate coefficients of K° and K % separately to zero. 
Prom the coefficient of K° we have 

2[(2» + <7o) 2 - (2/- 2* + l) 2 — 2m (2j - 2i + 1) -f m 2 ] iT i% _, 

This equation holds for all integral values of/, thus forming a system of equa- 
tions from which the quantities K may be determined. Now previously the 
quantities K were determined from the system 

(2j + gfK J =2^M i K i+i 

i 

for all integral values of j. The two systems must necessarily be equivalent. 
To show this, we have from the latter system 

(2» + gfK i a J _ i = 2a j _i^M f _ i K f , 

r 
hence 

2(2. + gtKfl,_ t = 2^K j ^a j _ i M f _ i = 2^K J %a j _ i M i _ r 

= 2X [coefficient of p-V + i in (JL + mA u] 

j' 

= 2-^ [coefficient of ^'- 2i+1 in D 2 u + 2mDu + m 2 u + ^ (u + *)] 
= 2 [(2i-2»+l) s +am(2y-2t + l)+|m 8 ]2r iS ._ i +^y^. a . . , 
Now this is the former system, and the equivalence is established. 
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Again from the coefficient of K*, 

2 [(2i+g f~(2j-2i+iy-2m {2j-2i+l)-im^E i {W) 3 _ i -^ m ^K i {k1i!) i _ j _ 1 
+ 2[— (2i+<, )+(2j+2i +l + 2g f+ 2m(2j—2i+l+2g )+im z ]K i (M) j + i 

+ X l(2i + g f - (2/ - 2% + If - 2m{2j— 2% + 1) — *m»] (#%>,._< 

— I m'S (&!/)&_,_! 

+ 2g y £(2i + g )K i a j _ i = 0. 

This is the system of equations from which g yl and the quantities (T&h!) must be 
found. 

(Jt?M) is at our disposal. We make it zero in order to simplify the equa- 
tions. This is equivalent to defining 2Ka as the coefficient of the principal 
periodic term in the value of z correctly as far as the cubes of the inclination. 
When K has been thus defined, its ratio to Delaunay's y is determinate to the 
order of the square of the inclination. Of this ratio the part dependent on m 
only has already been found. 

For a value of j other than zero, the only quantity with a coefficient of 
order zero is (& 2 &')j- The equation (j) therefore determines (Jt?h')j. Similarly 
the equation for which j = determines g yi . 

If the quantities (&%') be found to the fourth order, the corresponding terms 
of 2 are found correctly to the seventh order. This is the degree of accuracy to 
which we have gone in the first and second powers of the inclination; g yi can 
then be determined to the sixth order. This is one term beyond Delaunay's 
value. Before however g yi is obtained to the sixth order, it is necessary to find 
(M) to the sixth order. As this has not been done, g y % is here given to the fifth 
order only. 



Solving the equations, 






(»')* = 



2 7 i 2 8 

25 



2 6 

and 3 . . 51 . . 219 4 5751 B 

9+ = - -y m + ~W + ~W ~ ~W 
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If, after substituting the assumed value of z in the differential equation, we 
had equated to zero the coefficient of Z? j+1 ~ 8 , and then proceeded as before, we 
should have got ihe equations 

2[(2i + g Y-(2j+ 2i+ If- 2m(2j+2i+l)-^f]K t (M) i+j 
+ 2 [- (2i+g y+(2j-2i + 1— 2g f + 2m (2j-2i+l—2g )+ 5|?] K { (W)_ i+j 



2 



■3-1 



+ 2[(2i + g Y-(2j+ 2i + lf—2m(2j+2i+l)-^L](m) i a i+J 



+ 2g y ,X (2i + g ) K t a i+1 = 0. 



-^2(^) 4 «-i-,-i 



The quantities (WW), and also g y ,, may be found from this system of equations, 
which therefore serve to verify the values already found. Moreover, this system 
will give g y , to the sixth order, without its being necessary to calculate (hh) 
beyond the fifth order. 

The value of g y , so found is 

3 2 , 51 3 . 219 4 5751 , . 7115 6 

The arithmetic values of the quantities are as follows : 

^ = — .00806 633, 

(##)_!= —.01224 0, 

(jp#) x =—.00277 8, 

(##)_, = — .00008 0, 

(Wh% =—.00002 2, 

(##)_, = — .00000 1 . 

In this case, in the arithmetical work, the same equations are used as in the 
algebraical. 

The expression for the ratio of the motion of the Moon's node to the Moon's 
17 
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mean motion, so far as it consists of a series in m multiplied by the square of the 
inclination, has been computed by Delaunay to be 

. ( 3 - 27 „ 843 „ 7185 ,.1 
' \ 2 16 128 1024 ) 

where m! = — : — , and in this paper to be — 2x! — . E?. where 
1+ on 1+ffl 



1+ m' l r 1 + 

3 a , 51 3 . 219 4 

2 2* 2 7 2 10 2 12 



3 ,,51 3 . 219 4 5751 , , 7115 6 
9+ - — -S- m + "or m + "ST m — -^Io- m + -SET m » 



and where 

or numerically, 

^= — 0.00806 633, 

and the ratio of y 2 :^" 2 is 0.99974 65, hence 

— -Mx — JP= 0.00746 483. y 3 . 
1 + rn ' 

The numerical value of Delaunay's terms is 

0.00746 397 y 2 . 

The difference is therefore 0.00000 086 y 3 . 

Now Delaunay gives 0"7133 as the daily mean motion corresponding to 
these terms, and therefore the correction does not exceed 0"0001. This is the 
value of the residue of Delaunay's series that might be inferred from the series 
itself. The correction is small compared with the discrepancy that exists 
between theory and Observation. No large correction could have been expected 
from the sequence of Delaunay's terms. From Delaunay's complete expressions 
for the motion of the perigee and node (Comptes Rendus, LXXIV, p. 29), it 
appears that the series in m multiplied by y 2 is of more importance in the 
expression for motion of the perigee than for the node, whereas the series mul- 
tiplied by e 2 is of more importance for the motion of the node than for the 
perigee. 

Professor E. W. Brown has given (Amer. Jour., vol. XVII) a method of 
great analytical interest of calculating new parts of the motion of the perigee 
and node from the coefficients of the period terms already found ; for example, 
g y t can by Professor Brown's methods be deduced, by algebraical processes only, 
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without further reference to the differential equations, from the values of the 
coefficients a, K, (hh), (TeM), (h'W). The method is an extension of the process 
whereby Dr. Hill obtained the principal part of the motion of the perigee before 
calculating any periodic terms other than the variation. 
In the equation 

' (US + 2 2 )* 

we substitute 

3 

where g = g -\-g y!1 lP, and where, in the value of z, terms containing Sg in the index 
are omitted, as they do not affect what follows. Then equating to zero the 
coefficient of K 3 .£ )j+g , and employing the symbols M defined by 

— 00 

we have 

2 (ay + g ) g+ % + (ay + g«? (»'),- - as Jf,_< (^A/) 4 



= coefficient of JE*. f* + » in feV— 1 



«0 { («0 + *V).( S + **>) + A \ '■ 



Multiplying by Kj and summing for all values of j, the coefficient of (#W) S aris- 
ing from the second and third terms on the left-hand side is 

(2i + g f£ i -2XM j _ t K j , 
which is zero. 

Since — Kj is the coefficient of Etf-w-s i n the value of z y V— l/a , on the 

right-hand side we get half the coefficient of K* in the constant part of the value 

hi 1 

of -p-j -. \ - — s-y-5 , the presence of the factor one-half being due to 

«o 1 «o* + u y *s + s y M + z; \* 

the fact that the same terms arise on considering the coefficients of K 3 £~®~ g and 

K^ j+g . Hence g y , is given by the equation 

40 y2 2 (ay + g)*K P = const, coeff. of K* in -^ ^ - TTi . 

Were # 7 , determined in this manner, the approximations required to deter- 
mine the coefficients (&%') would be largely reduced. Moreover the equation 
that previously determined g yl would serve as a useful equation of verification. 
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But in the first place it may be noticed that this method will only give g yi alge- 
braically to the same order that the quantities (kk), (kk), (kk) have been calcu- 
lated, whereas the solution by continued approximation gives its value to one 
higher order. Possibly a similar objection would hold in a numerical applica- 
tion of the method. Again the arbitrary constants of integration arose from 
differential equations from which x, the mass of the Earth and Moon, had been 
eliminated. % is therefore a function of these constants. On p. 145, Amer. 
Jour., vol. I, Dr. Hill gives 

li=.x(n — n'f = n 2 <4 [1 — im? + . . . .] -3 , 

but in the higher approximation proceeded to in this paper some terms propor- 
tional to K* would have to be added to this value. These terms have not been 
calculated, and it would be unnecessary to calculate them unless it were to apply 
Prof. Brown's method of computing^,. It would be possible, of course, to 
replace the arbitrary condition (kk) = by the condition that (i shall contain no 
terms proportional to K 3 , but this is merely throwing a heavy piece of work 
upon an earlier portion of the paper. Moreover, when this has been done, or 
when the correction to p has been calculated (whichever alternative is adopted), 
there still remains the labor of multiplying several series of powers of £. The 
process of solving by continued approximation was indeed long, but it is 
extremely doubtful whether the practical application of Prof. Brown's method (to 
the motion of the node at any rate) would not be longer. It must be remem- 
bered, too, that the coefficients (Wk) would still remain to be calculated. 

Calculation of the Quantities (k 3 ). 

Substituting in the differential equation, and equating the coefficients of 
gij+i±3g to zer0) we obtajjj t^ equations 

2 t(' 3 9o + 2i) 2 — (2/— 2i + l) 2 - 2m (2/ - 2i + l) — f<] o,_, (¥) { 
t 

i 

+2 too + 2i) 8 — (2/ _ 2» + 1 + 2g f- 2m (2j— 2i + 1 + 2g ) - f m 2 ] K ( (kk)^ 

-fm 2 2-£(&n-i-i = 0, 
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and 

£ [(30 O + 2if — (2» + 2j + l) 2 — 2m (2» + 2j + 1) -f m 2 ] a i+i (^ 

i 

+2 [fo>+ 2*) 8 - (2/+ 2t+ 1- 2<7 ) 3 — 2m (2/+ 2»+ 1 - 2^ ) — f m 2 ] JT, (ft' #)i+-j 

-tm 2 2^TO-i-^-i=0. 

In the first of these equations (ft 3 },-, and in the second (ft 8 )-,-, are the only quan- 
tities that occur with coefficients of order zero. Hence the equations, as they 
stand, serve to determine (ft 3 ),- and (k s )_j respectively. 

In either equation/ can be given all integral values, positive and negative. 
Hence either system will determine all the quantities (7c?) ; the values of these 
quantities can then be verified by the other system of equations. 

In this way the following results have been obtained : 

m _ _ 9 2 , 117 3 _ 1607 4 



cn-.= -^ 4 



2 8 
81 
2 9 



The arithmetical results are given below. 

(W^i = — .06974 0, 

(7 £ 3 ) = + .00168 5, 

(& 3 )_ 2 = — .00159 9, 

(#)! = + .00001 4, 

(ft 3 ). 3 = — .00000 6. 

As before, the system of equations was first reduced to a linear system with 
numerical coefficients. The accuracy required in the coefficient was estimated 
from the approximate values of the quantities p as given by their algebraic 
values. 

This completes the calculation of the terms involving the first three orders 
of the inclination. 



